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Abstract 

Among the most fundamental observables of nucleon structure, electromagnetic form factors 
are a crucial benchmark for modern calculations describing the strong interaction dynamics of the 
nucleon's quark constituents; indeed, recent proton data have attracted intense theoretical interest. 
In this letter, wc report new measurements of the proton electromagnetic form factor ratio using 
the recoil polarization method, at momentum transfers = 5.2, 6.7, and 8.5 GeV^. By extending 
the range of for which is accurately determined by more than 50%, these measurements 
will provide significant constraints on models of nucleon structure in the non-perturbative regime. 
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The measurement of nucleon electromagnetic form factors, pioneered at Stanford in 
the 1950s, has again become the subject of intense investigation. Precise recoil polariza- 



tion experiments 



l| established conclusively that the proton electric form factor falls 



faster than the magnetic form factor G^.^ for momentum transfers > 1 GeV^, in dis- 
agreement with results obtained from cross section measurements j2"5|. Precise data to the 
highest possible are needed, for example, to test the onset of validity of perturbative 
QCD (pQCD) prediction. to~totic fan, facto, bchav.0. g contain Generalized 
Parton Distributions (GPDs) [7[, and to determine the nucleon's model- independent impact 
parameter-space charge and magnetization densities jsj. 

The effect of nucleon structure on elastic electron-nucleon scattering at a spacelike mo- 
mentum transfer = —Q^ < is described in the one-photon-exchange approximation 
by the helicity-conserving and helicity-flip form factors Fi{q^) (Dirac) and F2{q^) (Pauli), 
or alternatively the Sachs form factors, defined as the linear combinations Ge = Fi — TF2 
(electric) and Gm = F1+F2 (magnetic), where r = Q'^/AM'^ and M is the nucleon mass. Po- 
larization observables, such as the beam-target double-spin asymmetry 9| and polarization 



transfer 
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11| provide enhanced sensitivity to the electric form factor at large compared 



to cross section measurements, for which Gm becomes the dominant contribution. The po- 
larization of the recoil proton in the elastic scattering of longitudinally polarized electrons 
from unpolarized protons has longitudinal (Pi) and transverse (Pt) components with respect 
to the momentum transfer in the scattering plane llj. The ratio Pt/Pi is proportional to 



R 



rip 

^E 

r<p 



-Un-^ tan — 

^^Pi 2Mp 2 



(1) 



where /ip is the proton magnetic moment. Eg, is the beam energy, E'^ is the scattered e~ 
energy, 6^ is the e~ scattering angle and Mp is the proton mass. Because the extraction 
of from the ratio ([1]) is much less sensitive than the Rosenbluth method 12j to higher- 
order corrections beyond the standard radiative corrections [l^ , it is generally believed that 
polarization measurements provide the correct determination of G^ in the range where 
the two methods disagree. Previously neglected two-photon-exchange effects have been 



141 ]. and are a highly active area of theoretical 



shown to partially resolve the discrepancy 
and experimental investigation. 

The new measurements of G^/G^^ were carried out in experimental Hall C at Jefferson 



Lab. A continuous polarized electron beam was scattered from a 20 cm liquid hydrogen 
target, and elastically scattered electrons and protons were detected in coincidence. Typical 
beam currents ranged from 60-100 /xA. The beam helicity was reversed pseudorandomly at 
30 Hz. The beam polarization of typically 80-85% was monitored periodically using MoUer 
polarimetry [15i] . 

Scattered protons were detected in the Hall C High Momentum Spectrometer (HMS) [3], 
a superconducting magnetic spectrometer with three focusing quadrupole magnets followed 
by a 25° vertical bend dipole magnet, operated in a point-to-point tune. Charged particle 
trajectories at the focal plane were measured using drift chambers, and their momenta, 
scattering angles, and vertex coordinates were reconstructed using the transport matrix 
of the HMS. For this experiment, the HMS trigger was defined by a coincidence between 
the pair of scintillator planes just behind the drift chambers and an additional scintillator 
paddle placed at the exit of the dipole. The size of this new paddle matched the acceptance 
of elastically scattered protons. 

To measure the polarization of scattered protons, a double Focal Plane Polarimeter (FPP) 
was installed in the HMS detector hut, replacing the standard Cerenkov detector and rear 
scintillators. The FPP consists of two retractable 50 g cm~^ CH2 analyzer doors, each 
followed by a pair of large-acceptance drift chambers with an active area 164 x 132 cm^. 
The tracks of protons scattered in the analyzer material were reconstructed with an angular 
resolution of approximately 1 mrad. 

Scattered electrons were detected in a large-acceptance electromagnetic calorimeter (Big- 
Cal) positioned for each to cover a solid angle kinematically matched to the ~ 7 msr 
proton acceptance of the HMS, up to 143 msr at = 8.5 GeV^. BigCal was assembled 
from 1,744 lead-glass bars stacked in a rectangular array with a frontal area of 1.2 x 2.2 m^ 
and a thickness of approximately 15 radiation lengths. The trigger for BigCal was formed 
from analog sums of up to 64 channels, grouped with overlap to maximize the efficiency for 
electrons at high thresholds of nearly half the elastic e~ energy, used to suppress charged 
pions and low-energy backgrounds. The over-determined elastic ep kinematics allowed for 
continuous in situ calibration and gain matching. The primary trigger for the experiment 
was a time coincidence between BigCal and the HMS within a ±50 ns window. 

Elastic events were selected by applying cuts to enforce two-body reaction kinematics. 
The electron scattering angle 6e was predicted from the proton momentum pp and the 
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beam energy, and the azimuthal angle 0e was predicted from 0p assuming coplanarity of the 
electron and the proton. The predicted electron trajectory was projected from the interaction 
vertex to the surface of BigCal and compared to the measured shower coordinates. The small 
area of each cell relative to the transverse shower size resulted in coordinate resolution of 
5-10 mm, corresponding to an angular resolution of 1-3 mrad, which matched or exceeded 
the resolution of the predicted angles from elastic kinematics of the reconstructed proton. 

An elliptical cut {^x / Xmaxf + {^y / Umaxf' < 1 was applied to the horizontal and vertical 
coordinate differences (Ax, Ay), where {xmax^Umax) are the Q^-dependent, Sa cut widths 
used for the final analysis. An additional cut was applied to the proton angle-momentum 
correlation Pp—Pp{Op) which further suppressed the inelastic background. No cut was applied 
to the measured e~ energy, because the BigCal energy resolution was insufficient to provide 
additional separation between elastic and inelastic events. Figured] illustrates the separation 
of the elastic peak in the Pp — Pp{6p) spectrum using BigCal. 



o 
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FIG. 1. (Color online) Elastic event selection for = 8.5 GeV^. The momentum difference 
{Pp ~ Pp{Gp))/Po^ where po is the HMS central momentum, plotted for all events (black dashed), 
events passing the 3a elliptical cut (blue solid), and events failing the cut (green dotted). The 
estimated background (red dot-dashed) integrated over the final cut region (black vertical lines) is 
approximately 5.9%. 

The dominant background was hard-Bremsstrahlung-induced 7r° photoproduction, 7 -|- 
p — )■ 7r° + p, in the 2.3% radiation length cryotarget, with the proton detected in the HMS 
and one or two decay photons detected in BigCal. The kinematics of this reaction 
overlap with elastic ep scattering within experimental resolution for near-endpoint photons. 
The contribution of quasi-elastic Al(e, e'p) scattering from the cryocell windows was also 
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measured and found to be negligible after cuts. The total background including inelastic 
reactions and random coincidences was estimated as a function of Pp — Pp{Op), as shown in 
figure dl using a two-dimensional Gaussian extrapolation of the (Ax, Ay) distribution of 
the background into the cut region under the elastic peak. A Monte Carlo simulation of 
elastic ep scattering and vr^ photoproduction was performed as a check on the background 
estimation procedure. The two methods agreed at the 10% (relative) level for wide variations 
of the cuts. 

The angular distribution of protons scattered in the CH2 analyzers measures the polar- 
ization components at the focal plane. The polar and azimuthal scattering angles (f) of 
tracks in the FPP drift chambers were calculated relative to the incident track defined by 
the focal plane drift chambers. The measured angular distribution can be expressed in the 
general form, 

N^{p,^,^) = N^ [1 + (ci ± AyPf^n cos + 

{s^ T AyPin sin ^ + 

C2 cos(2(/}) + S2 sin(2v9) + . . . (2) 

where A''^'' is the number of incident protons in the ± beam helicity state, s{p,'d) is the 
fraction of protons of momentum p scattered by an angle Ay{p,i)) is the analyzing 
power of the P+CH2 reaction, and P/^^ and P^^^ are the transverse components of the 
proton polarization at the focal plane. Ci, Si, C2, S2, ■ ■ ■ are the Fourier coefficients of helicity- 
independent instrumental asymmetries, which are cancelled to first order by the helicity 
reversal. Figure [2] shows the measured helicity-dependent azimuthal asymmetry — = 
^^^+'^ — ^^^-^ ~ Ay [Pjj^^ cos — P^"^ sin yj] , where Ay9 is the bin width, summed over 
all p and the d range 0.5° < ^9 < 14° outside which Ay ^ 0. 

The extraction of Pt, Pi, and Pt/Pi from the measured asymmetry at the focal plane 
involves the precession of the proton polarization in the HMS magnetic field, governed by 



the Thomas-BMT equation 17|]. The rotation of longitudinal Pi into normal P/^^ allows 
the simultaneous measurement of Pt and Pi in the FPP, which is insensitive to longitudinal 
polarization. The unique spin transport matrix for each proton trajectory was calculated 
as a function of its angles, momentum, and vertex coordinates from a detailed model of the 
HMS using the differential-algebra based COSY software {isl. The polarization components 
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FIG. 2. (Color online) Helicity difference distribution /+-/_ for = 3.5 GeV^, 0.5° < -d < 14.0°. 
The data are fitted with /+ — /_ = a cos ip+bsin ip (solid curve), resulting in a = (0.16 it 1.19) x 10^^ 
and b = (-3.99 ± 0.12) x 10"^ (^Vn.d.f. = 0.67). 

at the target were then extracted by maximizing the Ukehhood function defined as: 

^ event 



C{PuPi) = n [l + he,Af{Sf,P, + S^iPi)cos^, 

i=l 

-he,Af{S^P, + S^P^ siny,, + A^l (3) 



where h is the beam polarization, S^^. are the spin transport matrix elements, = ±1 is 
the beam helicity, and Aq is the false asymmetry. 

The polarization of the residual inelastic background passing "elasticity" cuts was ob- 
tained from the rejected events using the same procedure, and used to correct the polariza- 
tion of elastic events. The acceptance-averaged fractional inelastic backgrounds for = 5.2, 
6.7, and 8.5 GeV^ were Ninei/{Ninei + Nei) = (1.12±0.16)%, (0.77±0.12)%, and (5.9±0.9)%, 
respectively. The resulting absolute corrections to R were AR = (8.4 ±1.5) x 10"'^, 
(7.5 ± 1.3) X 10-3, and (6.0 ± 1.3) x 10"^. 

Since the beam polarization and the P+CH2 analyzing power cancel in the ratio, there are 
few significant sources of systematic uncertainty in the results of this experiment. The most 
important contribution comes from the precession calculation. An excellent approximation 
to the full COSY calculation used for the final analysis is obtained from the product of simple 
rotations relative to the proton trajectory by angles x<t> iii the non-dispersive plane and xe 
in the dispersive plane. x</> = li^p4>bend and xe = I'^p^bend are proportional to the trajectory 
bend angles 4>bend and Ohend by a factor equal to the product of the proton's boost factor 7 
and anomalous magnetic moment Kp. The relevant matrix elements in this approximation 



are Syt = cosx<t>, ^yi = sinx^, S^^t = sinx^sinxe, and 5*^./ = - cos sin xe- These simple 
matrix elements were used to study the effects of systematic errors in the reconstructed 
kinematics. 

The error A(f)bend due to unknown misalignments of the quadrupoles relative to the HMS 
optical axis leads to an error 'jKpA(j)f,end on Pt/ Pi- This uncertainty was minimized through a 
dedicated study of the non- dispersive optics of the HMS following the method of [19] , setting 
a conservative upper limit of |A0| < 0.5 mrad, which is the single largest contribution to 
the systematic uncertainty in R. The contribution of uncertainties in the absolute central 
momentum of the HMS and the dispersive bend angle O^end is small by comparison. The 
extracted form factor ratio showed no statistically significant dependence on any of the 
variables involved in the precession calculation, providing a strong test of its quality. 

Uncertainties in E^,, E'^ and Of, make an even smaller contribution. Uncertainties in 
the scattering angles in the FPP were minimized by a software alignment procedure using 
"straight-through" data obtained with the CH2 doors open. False asymmetry coefficients 
obtained from Fourier analysis of the helicity sum distribution /+ + /_ were used to correct 
the small, second-order contributions to the extracted polarization components. The result- 
ing correction to R was small (|Ai?| < 0.007) and negative for each Q^. The correction 
procedure was verified using a Monte Carlo simulation. 



Ee, GeV 






R lb ARstat. i ARsyst. 


4.05 


60.3 


5.17 ± 0.123 


0.443 ±0.066 ±0.018 


5.71 


44.2 


6.70 ± 0.190 


0.327 ±0.105 ±0.022 


5.71 


69.0 


8.49 ± 0.167 


0.138 ±0.179 ±0.043 



TABLE I. Results for R = ^ipG^^/G^j^, with statistical and systematic uncertainties. E(, is the 
beam energy, Of. is the central electron scattering angle, {Q"^) is the acceptance-averaged , and 
AQ"^ is the r.m.s. acceptance. 



The results of the experiment are presented in table HI Standard radiative corrections to 
Pt/Pi were calculated using the code MASCARAD [l3|, found to be no greater than 0.13% 
(relative) for any of the three values, and were not applied. Figure |3] presents the new 
results with recent Rosenbluth and polarization data and selected theoretical predictions. 
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FIG. 3. (Color online) Upper panel: The proton form factor ratio fipG^/G^^ from this experi- 
ment (filled black triangles), with statistical error bars and systematic error band below the data. 
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231] (Cloet), 



2|] (Guidal), and 



211] (de Melo), [22] (Gross) 
251 ] (Belitsky). Lower panel: The same data and theory curves as the upper panel, expressed as 



Theoretical descriptions of nucleon form factors emphasize the importance of both baryon- 
meson and quark-gluon dynamics, with the former (latter) generally presumed to dominate 
in the low (high) energy limit. Recent Vector Meson Dominance (VMD) model fits by 



Lomon 



20| include p'(1450) and aj'(1420) mesons in addition to the usual p, u, and 0, anc. 



a "direct coupling" term enforcing pQCD-like behavior as — > oo. de Melo et al. 21] 
considered the non- valence components of the nucleon state in a light-front framework, using 
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Ansatze for the nucleoli Bethe-Salpeter amplitude and a microscopic version of the VMD 



model. Gross and Agbakpe 



22| modeled the nucleon as a bound state of three dressed 



valence constituent quarks in a covariant spectator theory. Cloet et al. [23| calculated a 
dressed-quark core contribution to the nucleon form factors in an approach based on Dyson- 
Schwinger equations (DSE) in QCD. The disagreement between this calculation and the 
data at lower is attributed to the omission of meson cloud effects. 

The Dirac and Pauli form factors are related to the vector (iJ) and tensor [E) GPDs 



through sum rules ^]. Guidal et al. [2J] fit a model of the valence quark GPDs based 
on Regge phenomenology to form factor data. In this model, the ratio F^jF^ constrains 
the X — )■ 1 behavior of -E, where x is the light-cone parton momentum fraction. When 
combined with the forward limit of H determined by parton distribution functions, the new 
information on E obtained from precise form factor data allowed an evaluation of Ji's sum 



rule |7| for the total angular momentum carried by quarks in the nucleon. 

The data do not yet satisfy the leading- twist, leading order pQCD "dimensional scal- 
ing" relation F| oc Ff/Q^ [6]. The modified scaling Q'^F^/Ff oc \n^{Q'^/K'^) obtained by 
considering the subleading twist components of the light-cone nucleon wavefunction 25| . 
with A = 300 MeV as shown in figure El describes the polarization data rather well. This 
"precocious scaling" of Ff/Ff is a necessary, but not sufficient condition for the validity 
of a pQCD description of nucleon form factors. Despite progress in calculations based on 
light cone QCD sum rules 26|, pQCD form factor predictions have not yet reached the level 



of accuracy of phenomenological models such as 20l-l22l |24| when applied to all four form 
factors (Ff'^), underscoring both the difficulty of predicting observables of hard exclusive 
reactions directly from QCD, and the strong guidance to theory provided by high quality 
data such as the results reported in this letter. 
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